Self-consistent non-Markovian theory of a quantum state evolution for quantum 

information processing 
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It is shown that the operator sum representation for non-Markovian dynamics and the Lindblad 
master equation in Markovian limit can be derived from a formal solution to quantum Liouville 
equation for a qubit system in the presence of decoherence processes self-consistently. Our formula- 
tion is the first principle theory based on projection-operator formalism to obtain an exact reduced 
density operator in time-convolutionless form starting from the quantum Liouville equation for a 
noisy quantum computer. The advantage of our approach is that it is general enough to describe a 
realistic quantum computer in the presence of decoherence provided details of the Hamiltonians are 
known. 
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Dynamics of a quantum system coupled to an environ- 
ment has been studied extensively for potential applica- 
tions to quantum computing and quantum information 
processing recently The key element of the stud- 

ies is the reduced-density-operator which is a solution 
to quantum Liouville equation (QLE). The QLE would 
involve Hamiltonians for systems representing qubits, 
reservoir, and mutual interaction between the system and 
reservoir that causes decoherence || . The presence of de- 
coherence would be the most important obstacle to the 
ideal operations of quantum gates or quantum channels 
||. To overcome this difficulty, the quantum error cor- 
recting codes [JtJ and the decoherence free subspaces ||] 
of multiple qubit systems have been suggested. For both 
quantum error correcting codes and decoherence free sub- 
spaces the knowledge of the reduced density operator of 
the qubit system is essential. 

Up to now the information about the reduced density 
operator is obtained from Lindblad master equation 
in Markovian approximation or an operator sum repre- 
sentation (OSR) [jl0| in the non-Markovian case which 
is also known as Kraus representation. Even though the 
OSR provides better information about the qubit system 
than the Markovian formalism, somewhat surprisingly, as 
pointed out by Bacon et al. pd[ |, the former is obtained 
in the language of gates, i.e., the unitary transforma- 
tion, rather than from the solution to the QLE itself in 
the Hamiltonian formulation. In other words, most of 
the proposals for quantum computers or quantum gates 
have assumed particular forms of the unitary transforma- 
tions beforehand. In our opinion, it would be desirable 
if there is a way to obtain the direct solution for the re- 
duced density operator from the QLE to model physical 
implementations of the quantum computers. 

The QLE is an integro-differential equation and, in 
general, it is nontrivial to obtain the solution of the form 



P 



(i) 



not necessarily a unitary transformation if one consid- 
ers an open system interacting with a reservoir in the 
presence of decoherence processes. Sometime ago we 
studied the time-convolutionlcss reduced density oper- 
ator formulation to model quantum devices jl2 13 and 
noisy quantum channels [T^ ] . In this theory the memory 
kernels of the Volterra-type integral equation are solved 
self-consistently using the superoperator formalism and 
it was shown that both non-Markovian decoherence pro- 
cess and renormalization of the memory effects can be 
incorporated. 

In this paper we formulate a general non-Markovian 
theory based on a QLE and show that the OSR for the 
non-Markovian case and the Lindblad master equation 
approach within the Markov approximation can be de- 
rived self-consistently. 

The Hamiltonian of the total system is assumed to be 



H t {t) = H s (t) + H b + H h 



(2) 



where H s (t) is the Hamiltonian of the system, Hi, the 
reservoir, and Hi nt the Hamiltonian for the interaction 
of the system with the reservoir. Note that the system 
Hamiltonian H s (t) may contain time-dependent external 
field terms to control the qubit system. The equation 
of motion for density operator p t of the total system is 
given by a QLE as 



j t Pt(t)=-i[H t (t),p t (t)} 



-iC t (t)p t (t), 



(3) 



where p is the reduced density operator and £ is the su- 
peroperator of linear mapping. The superoperator E is 



where Ct(t) — £s{t) + £b + £int is the Liouville oper- 
ator. The Liouville operators are in one-to-one corre- 
spondence with the Hamiltonians. Here we use a unit of 
H = 1. The reservoir is assumed to be in the thermal 
state. However, the assumption may be extended to any 
time-independent reservoir states which commutes with 
the reservoir Hamiltonian, i.e., CbPb = 0. In order to 
derive and to solve an equation for the system alone, we 
employ the projection-operators |l5|,|l(| that decompose 
the total system by eliminating the degrees of freedom for 
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the reservoir. Time-independent projection-operators V 
and Q are defined as 



VX=p b Tr b (X), Q=l-V, 



(4) 



for any dynamical variable X. Here Tr b indicates a par- 
tial trace over the quantum reservoir. The information of 
the system is contained in the reduced density operator 
p(t) given by 



P (t) = Tr b p t (t) 

= TY b P Pt (t). 



(5) 



After some mathematical manipulations, the time convo- 
lutionless equation of motion for Vptit) — p b p(t) is given 
byfM! 

^-Vp t {t) = -iVC t {t)Vpt{t) + iVC t {t) CAf(i) - 1) Pp t (t), 
at 

where 

AT x (t) = l + i f dr H(t, T)QC t (r)V G(t, r). (7) 
Jo 

The projected propagator Ti.(t,r) and the anti-time evo- 
lution operator G{t,r) of the total system are defined as 

H(t,T) = Texp^-iJ ds QA(s)Q| (8) 

and 

g(t,r) =T c exp|^ ds C t (s)\ , 

where T and T c denote the time ordering and the anti- 
time ordering operators respectively. The formal solution 
to Eq. (g) is given by Q 

Vp t (t)=U(t,0)Vp t (0) 

-if dsU{t,s)VC t {s){M{s)-l}Vp t {s), (9) 



with 

W(t) = 1 + z^ dsU 3 (t,s)Tr b {£ int Z(s)(l- Z(s)) -1 p b ) 

xTr, 0)K{t, s)Uo\t, 0) (1 - Z{t))- X p b } . (13) 

Here, we define 

Z(t) = l-M- 1 (t), 
U s (t,T) =Texp|-i J ds £ s (s) j , 
U (t,r) = exp{-i(t - T)C b }U s (t,T), 



(14) 
(15) 



(16) 
(17) 



and 



(Q)TZ(t,r) =T c exp jz^ ds Uq\s, 0)C mt U 



where the projected propagator U (t, t) of the system is 
defined by 



U(t,r) = Texpj-i J dsVC t (s)V^. 



(10) 



Eq. (^) can be put into time-convolutionless form by sub- 
stituting 



p t (s) = g(t,s)p t (t) 



(ii) 



and after some mathematical manipulations, we obtain 
the reduced density operator p(t), which is an exact so- 
lution to the QLE, given in the form of Eq. (0), 



p(t) = £(t)p(0) 

= W-\t)U s (t,0)p(0), 



(12) 



b(s,on, (is) 



where Uo(t,r) is the non-interacting time-evolution op- 
erator of the system and the reservoir and lZ(t, r) is the 
anti-time evolution operator of the total system in the 
interaction picture [|l3|,^4|| . 

It is straightforward to obtain the time-convolutionless 
equation of motion for a reduced density operator p(t). 
From Eq. (H), we get 



dt 



p{t)=-i£ g (t)p{t)+C(t)p(t), 



with 



C(t) = -iTr b {£ mt Z(t)(l - Z{t))- X p b } 



(19) 



(20) 



where C(t) is a generalized collision operator and we use 
an anzatz V 'CintP = which is equivalent to neglect 
renormalization of the unperturbed energy of the system 
& 

In the following, we first show that the time- 
convolutionless equation of motion ( |l9| ) becomes the 
Lindblad master equation in the Markov approximation. 
The lowest-order Born approximation, which is valid up 
to the order (Hint) 2 , is used subsequently. The effect of 
C(t) on p(t) up to the second-order expansion becomes 



C {2) (t)p(t) 







/ dr Tr 6 


Hint i 


h 





Hint(r - t),p b p(t) 



(21) 



where Hi n t(t) is the Heisenberg transformation of Hint 
defined by Uo(t)Hi n t- For the specific form of the inter- 
action Hamiltonian, we assume a Caldeira-Leggett-type 
model 17 ll| given by 



Hint = y^y a ®b a 



(22) 



where v a is the Hermitian operator acting on the system 
and b a = J2k(9akd\. + S^Ofe) is a fluctuating bosonic 



2 



quantum field whose unperturbed motion is governed by 
the harmonic oscillator Hamiltonian for the reservoir, 



H b {t) = y^^fcajafc- 



(23) 



The set of operators {v a } describes the various decoher- 
ence processes and sometimes they are denoted as the 
error generators. From Eqs. ( pi]) and (|22|), we obtain 

C (2) WpW = V I dr Xa p{T-t)[vp{T~t)p{t),v a ] 
a/3 Ja 

+ dr Xap(t - r)[v a , p(t)vp(r - t)} 

a/3 J ° 

(24) 



where 



Xapit) = Tr 6 6 a (i)6 /3 p fc = Tr b b a b/3{-t)p b . 



(25) 



The characteristic function Xapit) for the heat bath sat- 
isfies Xap(t) = X*p a (~t)- f n t ne Markovian limit, it be- 
comes 



(26) 



Then, we get 

C (2) (t)p(t)« \Y,1cct}{[v a m,vp} + [v a ,p{t)v fj }} (27) 

a/3 

where 7 Q/ 3 contains the information about the physi- 
cal decoherence parameters. It is now obvious that 
Eq. (|27|) is equivalent to the Lindblad term Cd described 
in Ref. ||, which takes into account the nonunitary, de- 
cohering dynamics. 

We now proceed to prove that the OSR or the Kraus 
representation can be derived from the formal solution 
given m Eqs. (|l|) and (|l|). The evolution superoperator 
£{t) becomes 



£^{t) = \l-i dsU s (t,s)Ti b C mt Z^(s)p b U^fas) 



The superoperator £^(t) satisfies the following con- 
ditions: (i) trace-preserving, (ii) Hermiticity-preserving, 
and (hi) complete positivity. As a result, there exists a 
corresponding OSR Jl(J. We will find the OSR for £ ^ (t) 
in Eq. (|2^) although any order of perturbation is appli- 
cable based on our formulation. Let {K a } be the set of 
Kraus operators for p(t) described in Eq. (|3C|), then 



p(t) = £ {2) (t)p{0) = K a {t)p{Q)Kt {t) 



(31) 



with the completeness relation, independent of the evolv- 
ing time t, 



Y,Kt{t)k a {t) = t. 



(32) 



In order to derive explicit expressions for the superoper- 
ator £( 2 \ we employ the interaction picture for the time 
evolution of the system state as 



p{t)=£^{t)p{0) 

= U- 1 (t,0)£^(t)p(0) 



(33) 



To derive the set of Kraus operators {K a }, we adopt a 
matrix representation for them. Then, 



£ {2) e -Vf'et 



lib 



with 



£ 



ah 



(e h £ {2) e ) 



(34) 



(35) 



where {e ab \e ab — \a)(b\} is an orthonormal basis set 
which spans the Hilbcrt-Schmidt space of reduced den- 
sity operators. The Kraus operator is expanded in this 
basis as 



xU s {t, 0) 



(28) 



with 



then, 



Z {1 \s) =-i drU a {s,T)L int U^{s,r), 



(29) 



within the Born approximation. Substituting Eqs. ( |22] ) 
and ph into Eq. (pi), Eq. (Til) becomes 



p{t) =U s (t,0)p(0) -U s (t, 0)J2 f ds f dr 

a/3 J * J ° 

x { Xaf3(r - s) [p(0)vp(T)$ a (s) - u Q ,(s)p(Q)u /3 (r)] 



S KgenrnK^ — ^ ] K a K a e ab . 

a aafc 

Comparing Eqs. ( |34| ) and (j37|), we obtain 

cab \ 1 an K bm* 

C/ nm / , "'a ""a 



(36) 



(37) 



(38) 



The conversion to £^ is straightforward since U s (t) is 
unitary. From Eqs. (pO) and (Bq), we get 



+X* a p(r " s) [v a (s)v p (T)p(0) - v (r)p(O)^a(s)] (30) 



C... 



b 

rim 



$an$< 



bin 



Ban St 



bm 



fianBbrn 



(39) 



where S nrl is a Kronccker delta, 
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Ban = T2 ds dr xtp(r- s)(a\v a (s)tip(T)\n), (40) 
a/3 J ° Jo 

Am,6m = V] / ds dr Xafi{r - s) 
„« Jo Jo 



If we define pit) 
p(t) becomes 



2Re/(i) - \f(t)\ 2 , then the OSR of 



p(t) = (l-p(t))lp(0)t+p(t)a z p(0)*z- 



(49) 



x(a\v a (s)\n)(b\vf}(T)\m)* . 

The set of Kraus operators are not unique and can 
be generated from a canonical set by an extended uni- 
tary matrix |l(],[l!|||,|o| . We will obtain the canonical set 
of Kraus operators. The superoperator can be re- 
garded as a positive and Hermitian matrix with (a, n) be- 
ing the row index and (b, m) the column index |2| . Then, 
there exists some unitary matrix U an ,a which diagonal- 



izes as 



cab _ tj j tt* 

L 'nm — / j u an,aU- a <~> om a . 



(42) 



Since all eigenvalues d a are positive, d a — y/d^y/d^, and 
Eq. ( fi"2| ) is in the form of Eq. (|3|) . One may choose k„™ 
as 



It is obvious that the Kraus operators ( ff8[ ) satisfy the 
(41) completeness relation (^). 

In summary, we have shown that the OSR for the non- 
Markovian case and the Lindblad master equation for 
the Markov case can be derived from the formal solution 
to the QLE for the qubit system in the presence of de- 
coherence processes self-consistently. Our formulation is 
the first principle theory starting from the exact solution 
to the QLE in time-convolutionless form and the ma- 
trix representation of the evolution superoperator. The 
advantage of our first principle theory is that it is gen- 
eral enough to model a realistic quantum computer in 
the presence of decoherence provided that details of the 
Hamiltonians for the system, reservoir, and the mutual 
interaction are known. 



da Uar. 



(43) 



All equivalent sets of Kraus operators for the given su- 
peroperator can be generated by "unitary remixing" of 
the canonical set with the eigenvalue vector d' extended 
by some arbitrary number of zeros as d' = (d, 0, ...,0) 



In addition to the derivation from the canonical set, 
when the superoperator is already in the form of Eq. (38), 
the Kraus operators can be obtained more explicitly. As 
an example, let us consider a simple dephasing channel 
for a single qubit system where Hamiltonian is given by 



where & = £kSk4 
(Ell), we obtain 



E 



(44) 



9 k a k . 



From Eqs. (§|), @, and 



X (t) = Tr b (b(t)bp b ) 

Ban ^an/(^)i 
Aan.bm = 2Re/(^)A a \b$an$bm 



(45) 
(46) 
(47) 



where A a is an eigenvalue of <r z , and fit) — 

Jo d s Jo ^ T X*( r — s )- Note that the eigenvectors |n) of 
a z was used for the basis to obt ain B and A. If we set 
«8 n = Ml ~ /(*)] and nf n = v/2Re/(i) - \f(t)\ 2 \ a 6 an , 
then the resulting Kraus operators are 



K, 



1 
1 



[1 -/(*)] 
V2Ref(t) - |/(i)| 2 



1 

-1 



(48) 
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